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ABSTRACT 

The detailed workings of the central engines of powerful quasars remain a mystery. 
This is primarily due to the fact that, at their cosmological distances, the inner re- 
gions of these quasars are spatially unresolvable. Reverberation mapping is now be- 
ginning to unlock the physics of the Broad Emission Line Region (BELR) in nearby, 
low-luminosity quasars, however it is still unknown whether this gas is dominated by 
virial motion, by outflows, or infall. The challenge is greater for more distant, powerful 
sources due to the very long response time of the BELR to changes in the continuum. 
We present a new technique for probing the kinematic properties of the BELR and 
accretion disk of high-z quasars using differential microlensing, and show how substan- 
tial information can be gained through a single observation of a strongly-lensed quasar 
using integral field spectroscopy. We apply this technique to GMOS IFU observations 
of the multiply-imaged quasar Q2237-I-0305, and find that the observed microlensing 
signature in the GUI] broad emission line favours gravitationally-dominated dynamics 
over an accelerating outflow. 
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1 INTRODUCTION 

Differential microlensing offers the promise of indirectly 
measuring the spatial size and kinematics of different emit- 
ting regions in quasars. In this paper, a new observation 
of differential microlensing between the continuum and dif- 
ferent velocity components of the GUI] broad line emission 
line in Q2237-I-0305 is presented. These data are then used 
to constrain the kinematics of the broad line gas. 

The core regions of quasars are unresolved at any wave- 
length, and so physical understanding must be based on 
theoretical models constrained by their observed spectral 
properties. The gas responsible for the broad emission lines 
presents a particular challenge as there are relatively few 
measurements capable of even beginning to probe its kine- 
matical structure. As such, there is still no agreement on 
even the simplest models of this structure. Beyond the obvi- 
ous width of the emission lines, which indicate large Dopplcr 
velocities of tens of thousands of km/sec, most of our cur- 
rent understanding has come from reverberation mapping 
and from polarisation observations of emission lines. The 
most recent reverberation signatures have been interpreted 
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as arising from outflowing, inflowing a nd virialised motions 
(|Dennev et aZ.ll2009l : ISentz et aLll2009l '). while the polarisa- 
tion studies have indi cated an outflowin g helical motion in 
the emission line gas (| Young et al\\2QQi ). Theoretical mod- 
els h ave favoured outflowing winds (eg iMurrav fc ChianS 
Il997l ). but altern ative points o f view are still under active 
consideration feg. lGaskelill2009l 'l. 

Differential microlensing provides an independent probe 
of the inner structure of quasars. In quasars that are subject 
to strong gravitational lensing, the BELR and potentially 
the accretion disk are resolved on the spatial scale of the 
fine magnification structure of the lens. Given a model for 
the lensing galaxy, observation of differential magnification 
between these components allows constraints to be placed 
on their sizes. 

The use of this technique was first discussed bv lNemirofll 

(!l988h . who specifically looked at the effect of a single low 
mass star on a range of kinematical models for the broad 
emission line region (BELR). A key result of this paper was 
that different parts of the emission lines were differentially 
magnified, depending on the kinematical model: in general 
the smaller spatial re gions showed the greatest variat ions, 
as would be expected. ISchneider fc Wambsgans"3 l)l99Cl l fur- 
ther considered the problem, noting that, in the case of a 
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macro-imaged quasar, the differences between the lines in 
different images can be used to test whether particular emit- 
ting regions are being significantly microlensed. In particu- 
lar, these authors noted that Keplerian motions in the BELR 
would be much easier to detect than infall (or presumably 
outflow) . A key improvement in the modelling compared to 
iNemirofj l|l988h was the use of microlensing magnification 
patterns to m odel the possib le statistical variations for each 
macro-image. ISchneider fc Wambsganss (.1990, ) also noted 
that differential microlensing would also effec t the redshift 
measured f rom a particular line. More recentlv lAbaias et al\ 
l|2007l ') and iLewis fc Ib"atal ( |2004l ) both considered different 
signatures which might be induced in the observed structure 
of the broad emission lines. 

There are at least four differential microlensing exper- 
iments which could allow the measurement of the physi- 
cal parameters of the quasar emission regions. Target-of- 
opportunity observations of a quasar c rossing a caustic pro - 
vide the cleanest imaging experiment (IWvithe et a/.ll2000l ). 
The physical interpretation of a caustic crossing event is 
straightforward. However in order to trigger the target-of- 
opportunity, regular monitoring is required; and for a rea- 
sonable annual probability of observing a caustic crossing, 
more than ten objects would need to be monitored. Recently, 
monitoring data has been used to fit the size of the region 
emitting t he quasar continuum using Baye sian Monte Carlo 
methods l|Kochanekl |2004| : lEigenbrod et a l. 2008b). These 
analyses rely on detailed modeling where a standard accre- 
tion disk is convolved with microlensing networks. A third 
and more specialized possibility arises when two macro- 
images straddle a caustic. In a surprising number of cases, 
the fluxes of the two images differ, while theory predicts 
them to be the same. Several variables in the modeling can 
affect the relative fluxes, but the size of the emi ssion re- 
gion is the dominant fact or l|Congdon et al. 2007; B ate et a/1 
l200'i1 ). lBate et al\ (|2007| ) have shown that measurements of 
the anomalous fluxes can be used to set limits on the size of 
the emission regions. 

The flnal method uses spectroscopic data to compare 
the shapes of emission lines and continuum spectra. If the 
broad emission line region has ordered kinematical motions, 
then differential magniflcation may change the shape of the 
emission lines of one image. In the case of a macro-imaged 
quasar, differences between the line spectra of the images 
can then provide a diagnostic for the kinematical motions 
of the broad emission lines. Such diffe rences have a l ready 
been observed in several sources (e.g. 'Keeto n et al. 20061: 



Eigenbrod et tt/.ll2008bl : ISluse et a l. 2008; Huts emekers et al\ 
2010t ). In this work we apply this last method to new inte- 



gral field spec troscopy (IFS) of t he multipy imaged quasar 
Q2237-f0305 l|Huchra et aLlll985l ). 

Q2237-I-03Q5 is an ideal candidate for differential mi- 
crolensing experiments. The close proximity of the lens- 
ing galaxy, at Zd — 0.0394 compared to the quasar red- 
shift of Zs — 1.695, yields a large projected Einstein Ra- 
dius (ER) of ~ 2 X 10" h^Q^ {M/MQ^^^cm in the source 
plane. This ER, which characterizes the size-scale of mag- 
nification fluctuations, is signiflcantly larger than its esti- 
mated continuum region size cm 
l|Witt. Mao fc Schechteilll995l ). and it is similar to the up- 
per size limit of its CIII] and Mgll BELR, determined from 
our previous observations of differential microlensing in this 



source l|Wavth et ffl/.ll2005l . hereafter W05). This indicates a 
high probability of differential microlensing between contin- 
uum and broad lines, and also a reasonable probability of 
differential microlensing within the BELR. 

Differential microlensing between the continuum and 
BELR in Q2237-I-0305 has been observed many times. Mi- 
crolensing wit hin the broad line i tself h as also been observed. 
Most notably, 'Eig enbrod et al\ l|2008h have conducted VLT 
monitoring of Q2237-I-0305 that spans several years, and 
present high quality slit spectroscopy of both continuum 
and BELR change, comparable to the data presented in this 
paper. However IFS observations can provide more reliable 
spectra as they allow more careful deblending of the small- 
separation lensed images and lensing galaxy. 

Interpretation of the differences between the spectra 
of the four images requires a model for both the accretion 
disk and the spatial and kinematical structure of the BELR. 
Broadly, larger flux changes are expected for regions closer 
to the accretion disk. The velocity structure, projected along 
the line-of-sight to the quasar should discriminate between 
models where the higher velocities are far from the disk (due 
to an accelerating wind) and those where a high velocity 
might be found close to the disk (due to Keplerian motion 
around the central black hole). It is still not known which 
(if either) of these two models best describes the kinematics 
of the BELR. 

The current work does not attempt to construct a fully 
consistent model of the BELR, nor include a full pho- 
toionisation calculation. Rather, a more limited question 
is addressed: near the accretion disk where the continuum 
is emitted, is the BELR gas moving relatively quickly or 
slowly? Answering this question will help resolve the domi- 
nant physical process in the inner BELR. 

In Section 2 of this paper we present the observational 
data, and the methods of data reduction and spectra extrac- 
tion. Section 3 describes the flux ratios for different emis- 
sion lines and at different wavelengths for the continuum. In 
Section 4 we present simple models to describe the BELR 
kinematics. Section 5 describes the computation of the mi- 
crolensing simulations which are convolved with the models. 
In Section 6, the results of comparing the model predictions 
with the data are presented. These are discussed in Section 
7, with the flnal conclusions presented in Section 8. 



2 OBSERVATIONAL DATA 

Q2237-f 0305 was observed on the 2 7th of June 2006 with 
the CMOS Integral Fie ld Unit (IFU) (|Allington-Smith et all 
l2002l : lHook et aZ.ll200j ) on the Gemini South telescope (pro- 
gram ID: GS-2006A-Q-14). 5x15 minute exposures were 
taken in one-slit mode using the B600 grating and with cen- 
tral wavelength 515 nm, giving a spectral rang e of 3975A to 
6575A and a resolution of 7? = 1688. At Q2237+0305's red- 
shift of 1.69, this covers the CIV and CIII] broad emission 
lines. Seeing conditions were excellent, with a PSF FWHM 
of 0'.'6. 

The CMOS IFU in one-slit mode gives a 5" x3'.'5 fleld of 
view, with an additional 5"xl'.'25 fleld offset by ~1 arcmin 
to measure the sky. Importantly, the field is 'filled', with 
the science field consisting of 500 0'.'2 diameter hexagonal 
lenslets with abutting edges. The filled field means that no 
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Figure 1. 2-D image created from collapsing ADR-corrected IFU 
data cube. Letters indicate the standard designation of quasar 
images in the Einstein Cross. 



interpolation between lenslets is required, which makes the 
determination of brightness ratios significantly more accu- 
rate than in the case of non-filled IFUs. 



2.1 Data Reduction 

The data was reduced using the Gemini IRAF package vl.8, 
under IRAF v2.12. The standard CMOS IFU reduction 
recipe was followed for overscan subtraction and trimming, 
bias subtraction, dark subtraction, flatfielding, arc calibra- 
tion, spectrum extraction, sky subtraction, and fiux calibra- 
tion. We used calibration files taken as part of the standard 
Gemini Facility Calibration Unit (GCAL) suite. 

Cosmic ray removal was performed separately on the 
reduced images before spectrum extraction, using CRREJ 
algorithm via IMCOMBINE. CRREJ performs iterative re- 
jection of unusually high pixels, accounting for the noise 
properties of the image. In W05 we found that this method 
provided far superior cosmic ray rejection to that provided 
by GSCRREJ, and demonstrated the advantage of taking 
multiple GMOS IFU exposures at the same offset position, 
as opposed to using an extensive dither pattern. We con- 
verted the hexagonal image arrays into rectangular arrays 
with (/.'l square pixels using GFCUBE. 

Atmospheric differential refraction (ADR) was found 
to cause a 0'.'8 shift over the wavelength range. This was 
corrected by first computing the shifts relative to the 



image position at the central wavelength of CIII] using 
cross-correlation, then fitting a physical model of ADR 
(jSchubert fc W alterscheidI l2000l ) to these shifts. The cube 
was corrected according to the best-fit ADR model. 

Figure [1] shows the ADR-corrected data cube collapsed 
into a 2-D image. The four quasar images and the core of 
the lensing galaxy are well-resolved. 



2.2 PSF Construction and Extraction of Image 
Spectra 

The excellent seeing allows us to extract the spectum of each 
quasar image with high precision. Nonetheless, there is some 
blending between adjacent quasar images and with the core 
of the lensing galaxy, and so deblending is required. 

To deblend, we perform least-squares fitting of the spa- 
tial frame at each wavelength element. With positions de- 
termined from archival HST imaging, we fit a reconstructed 
Point Spread Function (PSF; see below) to all quasar images 
and the galaxy core using minimization. For each quasar 
image, we perform PSF-subtraction of all other quasar im- 
ages and the galaxy core using the best-fit amplitudes. We 
then determine the fiux for each image, at each wavelength 
element, within an integration aperture, normalized by the 
fraction of the model PSF within this aperture. This aper- 
ture is a O'.'S radius circle centred on the relevant quasar 
image or galaxy, excluding regions within O'.'S of any other 
quasar image or the galaxy core. 

The PSFs used in these fits are determined using an 
iterative process. As the PSF is wavelength-dependent, we 
calculate ten different PSFs across the wavelength range, 
each from a ~100A wide stacked frame. The PSF extraction 
proceeds as follows: 

a. Take the azimuthal average for each quasar image 
within l'.'6 circles, excluding regions within O'.'S of any other 
quasar image or of the galaxy core. This yields four 1-D PSF 
profiles. We splice these together to construct a single PSF 
profile and convert this into a azimuthally symmetric 2-D 
PSF. 

b. Perform least-squares fitting to obtain a first estimate 
of quasar image and galaxy core amplitudes. 

c. Create four frames; in each one, three of the quasar 
images and the galaxy core are subtracted according to the 
best-fit normalizations from step h, leaving a single quasar 
image remaining. 

d. Use these four frames to splice together a 2-D es- 
timate of the PSF using the same regions defined in step 
a. 

e. Repeat steps c through d using the 2-D PSF. 

/. Repeat until the PSF converges, with RMS differences 
between iterations of < 0.1%. 

Using this process, the PSF converged within 5 itera- 
tions for all ten lOOA wavelength bins, and was stable for 
further iterations. 

Figure [2] illustrates the regions used for construction 
of the PSF and spectral extraction (left), and the residuals 
after final subtraction for the spectral region of the CIII] 
line. Figure [3] shows the extracted spectra of the four quasar 
images. 

The use of the PSF to model the core of the lensing 
galaxy introduces some error. If we assume a d e Vaucouleurs' 
effective radius of 3'.'l (|Schmidt et al\ll99^ . we find that 
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Figure 2. Left: Each spectrum is extracted from a O'.'S aperture centred on the relevant quasar image or the lensing galaxy core (solid 
circles), but excluding regions overlapping any other O'.'S aperture. The model PSF is constructed with contributions from each quasar 
image, using 1''6 apertures (dashed circles), and again excluding regions within 0''8 of other images or the galaxy core (see Sect. [T2l for 
details). Right: the residual map around the CIII] line after the best- fit subtraction of the constructed PSF from each quasar image and 
galaxy core. The median of the absolute noise-weighted residuals is ~ 0.2(T, with a range ol ± 2a 
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Figure 3. Extracted spectra for the four quasar images and the lensing galaxy in Q2237+0305: black: image A, red: image B, green: 
image C, blue: image D, orange: lensing galaxy. Dashed lines mark positions of CIV and CIII] broad emission lines. Solid vertical lines 
indicate positions of galaxy absorptions lines, from left to right: Mgb, Fell, and NaD 
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the difference in galaxy flux between the de Vaucouleur and 
PSF models in the integration aperture of images A and B 
is less than 1% of the quasar image flux in this region. As 
the galaxy spectrum is reasonably smooth in the region of 
the CIII] line, any residual galaxy flux is largely removed in 
the continuum subtraction (see Section [3. 3p . 



3 FLUX AND SPECTRAL RATIOS 

In a multiply-imaged quasar, differential microlensing is 
characterized by a variation in the magnification levels of 
different spectral components in a single lens image, com- 
pared to the "true", unmicrolensed quasar spectrum. This 
effect is most easily observed in the ratio of lens images. 
In the absense of microlensing, we expect emission line and 
continuum flux ratios to be equivalent to each other and 
to agree with the macro-model flux ratios, and we also ex- 
pect flat, featureless spectral ratios. Gravitational lensing of 
a point source is achromatic, and so the spectral differences 
resulting from differential microlensing can be identifled in- 
dependently of the lens model. However it is important to 
consider the other processes that can also lead to spectral 
differences between lens images. 



3.1 Differential Extinction and Intrinsic 
Variability 

It is possible that extinction properties differ signiflcantly 
between the lines of sight of each quasar image; both within 
the lensing galaxy and in intervening absorbing systems. 
However in this analysis we are concerned with the rela- 
tive strengths of the continuum and BELR within a small 
wavelength range (the width of the emission line), and 
so differential extinction will not be a signiflcant effect. 
Nonet heless, we correct fo r extin ction using the Ay val- 
ues of lAgol. Jones fc BlaesI (|2000l ). which uses the average 
Milky Way extinction curve to obtain Av for images A 
of 0.88±0.21, B of 0.84±0.20, C of 1.30±0.31, and D of 
1.15±0.27. 

Intrinsic variation within the quasar coupled with time 
delays between lens images can result in time-varying spec- 
tral differences, possibly duplicating the effect of differential 
microlensing if the continuum undergoes signiflcant varia- 
tion over a timescale smaller than that of the time delay. In 
the case of Q2237-f0305 the time delay pre dicted by mod- 
eling is less than a day (|Kent fc Falcdll988l ). and Chandra 
observatio ns measure a de lay between images A and B of 
~ 3 hours l|Dai et oZ.ll2"00^ . This is sufficiently smaller than 
the estimated accretion disk light-crossing time of ~ 6 days 
to eliminate the possibility of significant spectral differences 
due to intrinsic variation. 
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Figure 4. Ratios of spectra for images B, C, and D to image A. 
Residual features at locations of broad emission lines CIV and 
cm] (dashed lines) reveal differential microlensing in which the 
continuum in image A is magnified relative to the BELR. 



Table 1. Flux ratios for CIII] and CIV broad emission lines and 
for the continuum at the central wavelengths o f these lines. Mid 
IR flux ratios from lAeol. Jones fc BlaesI I2OO0II, radio flux ratios 
from lFalco et al\ lll996h . Hfl from lMetcalf et al\ ||2004|) . 



Region 


B/A 


C/A 


D/A 


CIII] line 


0.569 ±0.022 


0.643 ± 0.025 


0.630 ±0.026 


CIII] cont. 


0.431 ±0.026 


0.443 ±0.027 


0.281 ± 0.018 


CIV line 


0.473 ± 0.00; 020 


0.454 ±0.020 


0.422 ± 0.021 


CIV cont. 


0.375 ± 0.034 


0.347 ±0.033 


0.203 ± 0.025 


Mid IR 


1.11 ±0.11 


0.59 ±0.09 


1.00 ±0.10 


8 GHz 


1.08 ±0.27 


0.55 ±0.21 


0.77 ±0.23 


Hp line 


0.376 ±0.007 


0.387 ±0.007 


0.461 ± 0.004 



3.2 Flux Ratios 

Table [1] gives the flux ratios in the CIV and CIII] emission 
lines and in the continuum. Line ratios use the integrated, 
continuum-subtracted line fluxes, while continuum ratios use 
the the continuum level at the central wavelength of each 
line, determined by a linear interpolation from either side of 
each line. From these ratios it is apparent that image A is 
magnifled relative to the other images. 



We also see differential microlensing between the BELR 
and the continuum, with the continuum magnified to a 
greater extent than the BELR. This is expected: given any 
level of differential microlensing, the smaller emission re- 
gion should exhibit the most extreme deviation from the 
macro-model magnification. Using a method similar to that 

used in W05, we find a size limits for the CIII] BELR of 

1/2 

< O.OGG/iyg pc and for the continuum emission region of 

< 0.059^7g^ pc, both at the 95% confidence level. Note that, 
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while the limit on the BELR size is similar to that found 
from our 2002 data, the continuum limit is less tight than the 

1/2 

limit of < Omhj'ff pc reported in W05. This is because the 
earlier microlensing event resulted in stronger magnification 
of the continuum — possible only for smaller accretion disk 
sizes. Our continuum size limit is also consistent with, but 
less ti ght t han, the limits found by Pointdexter fc Kochanelj 
l|201(][l and lEigenbrod et al\ (|2008bf K 

We omit the CIV line from this analysis due to the 
difficulty in performing an accurate continuum subtraction 
at the far blue end of the spectrum. 



3.3 Spectral Ratios 

Figure |4] shows the spectral ratios of each quasar image rel- 
ative to image A. The presence of features around 1549A 
and 1909A illustrates the differential microlensing between 
BELR and continuum. The narrowness of these features sug- 
gest that there may be differential microlensing across the 
velocity structure of the line. In the case of a static mag- 
nification difference between line and continuum we would 
expect a broad feature covering the extent of the emission 
line. We also note that the continuum spectrum of image A 
is bluer than that of the other images, indicating possible 
differential microlensing across the accretion disk. 

To investigate differential microlensing within the broad 
line region, we look at the continuum-subtracted emission 
line ratios. The continuum spectrum at each line is approx- 
imated as the linear fit to a 120A-wide region centered on 
the line, masked in the central 80A-wide region. Figure [S] 
shows the continuum-subtracted CIV and CIII] lines for im- 
ages A and B (upper plots) and the continuum-subtracted 
line ratios (lower plots), which are equivalent to emission 
line magnification ratios. 

Errors in the spectral line ratios are derived from un- 
certainties in the PSF fits (see Sect. 12. 2p . which take into 
account errors in subtraction of other images and of the 
lensing galaxy. We include the error in continuum subtrac- 
tion, determined conservatively from the error in the linear 
fit around each line — both from the fit itself and from 
the variation in the continuum level from fitting to differ- 
ent spectral regions. For the simple analysis to follow, we 
don't directly calculate the error due to contamination by 
any iron lines beneath the broad emission lines. Instead we 
assume that this uncertainty is small compared to the un- 
certainty in the continuum. We have no direct means to test 
this assumption. 

Both the CIV and the CIII] line ratios reveal a sym- 
metric trend in the magnification ratio about the line cen- 
tre. This trend is highly significant for the CIII] line. While 
the trend is apparent for the CIV line, the difficulty in con- 
tinuum subtraction for this line (see Sect. 13.211 dramatically 
increases the uncertainty in the spectral line ratios, and so 
we restrict further analysis of this differential microlensing 
feature to the CIII] line. 

At low velocities near the line centre, the CIII] broad 
line in image A is magnified compared to the macro-model 
expectation: B/A ~ 0.65 compared to the expected B /A ^ 
1.0. As line velocity increases this magnification ratio drops 
further, towards the continuum ratio of B/A ~ 0.4. The 
observed trend suggests that the lower velocity broad line 
emission samples a region of the caustic structure closer to 



the average macro-model magnification level for image A 
than both the higher velocity gas and the continuum. This 
in turn suggests that size of the emission region approaches 
that of the continuum region with increasing gas velocity. 
Note that this feature cannot be explained by differential 
microlensing between the narrow line region versus the BLR 
as the width of the feature exceeds 20 A, or ±1500 km/s. In 
the next section we analyse this microlensing signature in 
the context of simple models of BELR kinematics. 



4 MODELS OF THE BELR 

We test simple models for two of the physical processes that 
are considered to be the primary candidates as the domi- 
nant drivers of BELR kinematics: an outflow driven by ra- 
diation pressure and orbital motion driven by gravity. The 
aim of these models is to simulate the expected wavelength- 
dependent surface brightness distributions for these models, 
and then, by convolving these with simulated magnification 
maps for Q2237-I-0305, to determine the expected magni- 
fication ratios as a function of wavelength, which we will 
compare to the observed microlensing signature. 



4.1 The Outflow Model 

For our simple outflow model we assume a radial wi nd in 
twin conical shells. Drawing on the wind model of lElvid 
(|2000l) . we consider optically thick clouds driven outwards 
by radiation pressure, with each cloud sufficiently small 
cross-sectional area to allow us to ignore self-shielding. 
Under these a s sumpt ions we may use the derivation of 
ICapriotti et al\ (|l980l ). which gives the acceleration of a 
cloud due to radiation pressure as: 



AcL 



(1) 



A-Kcr^Mc 

where Ac is the cross-sectional area of the cloud, Mc is its 
mass, L is the total ionising luminosity, and r is the distance 
from the ionising source. Velocity as a function of radius is 
therefore: 

AcL 



vdv ■ 



-dr 



(2) 



47r cr 2 M, 

So for a starting velocity t;o at an initial radius ro we have: 



v{r) 



ro 



(3) 



where K — AcL/2ncMc. And the terminal velocity, at r — >■ 
oo, is is VT = y^{K + rovD/ro- 

At a given projected distance from the centre, x, we 

know vios ~ v\J 1 — fr- So from equation [S] the radius at 
which a given vios occurs along a given line of sight satisfies 
the cubic: 



/(«) = a J + (- 







(4) 



where a = Ijr. This function yields a single physical solu- 
tion for the distance from the ionising source of clouds with 
projected velocity Vioa at a given projected distance from 
the centre, x. 
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Figure 5. Upper: continuum-subtracted CIV and CIII] emission lines for images A {solid) and B (dotted). Lower: B/A ratios using 
continuum-subracted emission lines. Solid lines show regions with significant line flux, while dashed lines show regions heavily dominated 
by noise. Dotted line shows the B/A continuum flux ratio. 



We wish to determine the expected surface brightness distri- 
bution produced by such a wind as a function of wavelength. 
If we assume that broad hne clouds re-emit a constant frac- 
tion of the incident ionising flux in the emission line, then 
the surface brightness of an emission line at AA from its rest 
wavelength and at projected distance x from the centre is: 



dSx^Sz 
tJ-x,A\ oc ^ — 



(5) 



where d is the number density of broad-line clouds at r, 5x 
is the size of the spatial bin, and 5z is the distance along the 
line of sight over which the emission line is shifted into our 
spectral resolution bin, assumed to be small compared to r. 
For a spherical outflow, if we assume cloud conservation (i.e. 
their lifetime is long compared to the flow timescale), then 
d satisfies: d oc l/(«r^). The size of projected spatial bins, 
5x, is constant, and so: 



5z 

fJ.x,AX OC j 



(6) 



So combined with equation[3]we have surface brightness 
at X and AA as a function of the distance from the source 
at which that emission is produced. For the purpose of the 
model it is useful to parameterize K in terms of ro, vq, and 

Sz 



A»a:,AA OC 



(7) 



where r and 5z are calculated using equation 2] 

We include both axial and equatorial obscuration with 
varying opening angles to simulate constrained winds or ob- 
scuration by jet and disk/torus. The final model of a three- 
dimensional axially symmetric radial wind is parameterized 
by: launching radius, terminal velocity (which encompasses 



cloud properties and continuum luminosity), orientation of 
the axis to the line of sight, and the opening angle of axial 
and equatorial black-out regions. 

Figure [6] shows an example set of surface brightness 
distributions for the wind model. As can be seen, significant 
asymmetries arise for axis orientations not perpendicular to 
the line of sight. Any obscuring disk or jet quickly blacks 
out the entire receding wind as velocity increases, and this 
obscuration increases with smaller orientation angles. Inter- 
estingly, the size scale of the surface brightness distribution 
does not change dramatically with gas velocity, although 
the red- and blue-shifted high velocity gas is offset in the 
plane of the sky for low orientation angles. Such offsets are 
likely to lead to strong differential microlensing, and hence 
asymetries between red- and blue-shifted line ratios. 



4.2 The Orbital Model 

To study the expected microlensing signature of an orbital 
BELR, we assume a purely Keplerian model consisting of 
randomly-oriented circular orbits. In this case the analog of 
equation |4l 



2 

-C / \ „2 3 2 ^los r\ 

f{a) = x a -a - "^=0 



(8) 



and as previously x is the 



where a — -iy, K' = 
projected distance from the centre and vios is the line-of- 
sight velocity that we're interested in. The solution to this 
cubic gives us the radius at which vios occurs for the current 

X. 

We follow the same approach as for the wind model, ex- 
cept that the cloud number density can't be constrained by 
simple conservation arguments, and so we need to choose a 
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Figure 6. Outflow model surface brightness maps for a selec- 
tion of the explored parameter space. The depicted models have 
launching radius ro = 1 X 10^^ cm and a range of disk/jet open- 
ing angles and orientations of the axis to the line of sight. For 
each parameter set, a small selection of the calculated velocity 
slices are shown, with negative velocities indicating the receding 
outflow. The surface brightness peak of the receding outflow is 
increasingly displaced relative to the approaching outflow as ve- 
locity increases. An obscuring disk and jet blacks out a larger 
fraction of the receding wind as orientation angle decreases and 
velocity increases, completely blacking out even lower velocity gas 
at angles below 60°to 70°. Each frame size is 3 X 10^^ cm. 

density gradient, d oc r" . We choose a value of a = —1, in the 
mid-r ange of the mid-range of likely values (e.g. iGoad et al\ 
1 19931 ). 

The surface brightness distribution of a purely Keple- 
rian BELR is then: 

^J.x,A\ <x K'Szr"'^ (9) 

The final model of a spherical BELR comprised of ran- 
dom circular Keplerian orbits is parameterized by: black hole 
mass, density gradient, orientation of the axis to the line of 
sight, and the opening angle of axial and equatorial black- 
out regions. 

Figure [7| shows an example set of surface brightness dis- 
tributions for the orbital model. This model produces sym- 
metric surface brightness distributions for the red- and blue- 
shifted gas. 

4.3 Parameter Range 

For both the outflow and the orbital model we include obscu- 
ration caused by an equatorial disk and an axial jet. We ex- 
plore the following disk/jet opening angles for these: 0°/0°, 



Figure 7. Orbital model surface brightness maps for a selection 
of the explored parameter space. The depicted models have black 
hole mass ttibh = 1 x 10^ Mq and a range of disk/jet opening 
angles and orientations of the axis to the line of sight. For each 
parameter set, a small selection of the calculated velocity slices 
are shown. Receding and approaching gas produces symmetrical 
surface brightness distributions. Each frame size is 3 X lO^'^ cm. 

0710°, 1070° , 10710°, 40°/40°, T/0°, T/10°. "T" indi- 
cates a transparent accretion disk; in all other cases the disk 
is opaque, meaning that only the near side of the BELR is 
visible. 

For both the kinematic models and all opening angle 
combinations we explore the following orientations: 0° , 15° , 
30°, 45°, 60° 70°, 80°, 90°; 

For the outflow model we simulate wind launching radii 
ro of 1 X lO^'^ cm, 3 x lO^'^ cm, 1 x 10^*^ cm, 3 x 10^® cm, and 
6 x 10^^ cm. We assume an initial velocity of uo = and a 
terminal velocity for the wind of Vt ~ 0.1c 

For the orbital model we simulate black hole masses 

rUBH of 1 X 10* Mq, 3 X 10* Mq, 1 X 10^ Mq, 1.5 X lO'' Mq, 

2 X 10^ Mq, 3 X lO'' Mq, 4 X 10^ Mq, and 6 x 10^ Mq. We 
assume a density gradient of a = — 1. 

For each model and parameter combination we calculate 
surface brightness maps for 5A-wide wavelength bins around 
the CIII] central wavelength between 1809A and 2009A. 



5 MICROLENSING SIMULATIONS 

To determine the microlensing signatures predicted by these 
BELR models, we perform a Monte Carlo simulation in 
which we convolve the model surface brightness maps with 
simulated magnification maps to produce a large set of ex- 
pected line ratio spectra. 
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The magnification maps ar e constructed using an in- 
verse ray-shooting tectmique (e.g. Kavser. Refsdal fc StabeUl 



ll986l : IWambsganss. Paczvhski fc Katj|l990al ). The key pa- 
rameters for such simuiations are the convergence utot 
and the shear 7 of the lens at the image positions. The 
Q2237-I-0305 lens model is provided in Tab le H (Trott 
fc W ayth, priv. comm., based on models of iTrott et al\ 
|2010|) . We assume a smooth matter component in the 
lens of 0%, and so the convergence is provided solely by 
a clumpy stellar component. This is reasonable as the 
lensed images in Q2237-I-0305 lie very close to the bulge 
of the lensing galaxy, where the stellar component is as- 
sumed to dominate ( Kent fc Fa lco 198^: [Schneider et al\ 
1 19881 : ISchmidt. Webster fc Le^.1998 : .Trott et 0^.1120101) . A 
fixed microlens mass of IMq was used. The Einstein Radius 
changes as the square root of this mass, although clustering 
of caustics is relatively unaffected. 

We generated magnification maps for images A and B. 
These maps have 2048 x 2048 pixels, with 5 different pixel 
scales for each. The scale size corresponds to the scale of 
the BELR model - which increases with ro and Mbh - 
and is chosen to allow sufficient independent positions for 
the Monte Carlo simulation (see below). Table [3] lists the 
pixel scale used for the different black hole masses and wind 
launching radii. 

For each BELR model parameter set, we generated sur- 
face brightness maps for 38 wavelength bins spanning 2OOA 
centred on the emission line central wavelength. At each it- 
eration of the Monte Carlo simulation we choose two source 
positions, one on the image A magnification maps and one 
on the image B magnification maps. We then convolve the 
BELR surface brightness maps for all 38 wavelength bins 
with the magnification maps, centred on these two posi- 
tions. The resulting pair of magnification spectra yields our 
magnificiation ratio spectra for the BELR model, which is 
equivalent to a continuum-subtracted line fiux ratio. 

The two magnification map positions are constrained 
by the observed continuum flux ratio. We first randomly se- 
lect a source position on the B magnification maps. Source 
positions on the A magnification map are then randomly se- 
lected with the constraint that the B/A magnification ratio 
for a point source be within 10% of the observed continuum 
flux ratio B/A = 0.39. 

The above assumes that the dominant continuum emis- 
sion region is smaller than the spatial resolution element 
at each pixel scale. For the smallest pixel scales (map 
sizes 2.4r;o and 7.2r]o, corresponding to pixel scales of 2.1- 
6.3x10^/1^/^ pc), this gives a continuum size smaller than 
the upper limit found in Section 13.21 However, these sizes 
are appropriate to the scale of the BELR models consid- 
ered. They assume the dominant optical continuum emis- 
sion arises from a region at a factor of ~ 3 or more smaller 
than the smallest scale of broad-line emission. 

We generated 1215 line ratio spectra per wind model, 
4860 line ratio spectra for the Mbh = 1.5 x 10^ Mq, 2 x 
10^ Mq and 4 x 10^ Mq orbital models, and 1215 line ratio 
spectra for all other orbital models. These line ratio spectra 
could then be compared directly with the observed B/A line 
ratio spectrum. 

Figures |S] and H] shows a random selection of simulated 
magnification ratio spectra for the outfiow model and the or- 
bital model respectively, for a single parameter combination. 



Table 2. 


Lensing 


parameters 


Image 




7 A^tot 


A 


0.413 


0.382 5.034 


B 


0.410 


0.384 4.984 



Lensing parameters for images A and B in Q 2237+0305 (T rott 
& Wayth, priv. comm., based on models of lTrott ei "aPbOlOl '). 



Table 3. BELR models and magnification maps 



Map size 



orbital models 
(Mf,h in Mq) 



wind models 
(ro in cm) 



2.4»7o X 2.4»7o 
7.2r;o X 7.2?7o 
24r)o X 24?7o 
72»)o X 72»7o 
144»)o X 144»)o 



1 X 10** 
3 X 10** 
1 X 10^, 1.5 X 10'' 
2 X 10^, 3 X lO'', 4 X 10^ 
6 X 10^ 



1 X 10^^ 
3 X lO^^^ 
1 X lO^'^ 
3 X lO^'^ 
6 X lO^'' 



Magnification map size, and the corresponding orbital and wind 
models, -qo is the Einstein Radius projected on to the source 
plane: 1.79 X 10^'' h~^^^ {M /MccYf^ cm for Q2237+0305. 



For outfiow model the parameters are: ro = 1 x 10^® cm, ori- 
entation 70°, disk/jet opening angle 0°/0°. For the orbital 
model the parameters are: Mbh = 3 x l(f M©, orientation 
0°, disk/jet opening angle 0°/0°. 



6 MEASUREMENT OF BELR KINEMATICS 

As illustrated in Figures [8] and |9l the orbital model tends to 
reproduce the microlensing signature observed in the data 
more easily than the wind model. To quantify this compari- 
son, we adopt a simple parameterization for the observed mi- 
crolensing signature, comparing the amplitude and width of 
the observed hump. Taking the CUT] continuum-subtracted 
line ratios between images A and B, we compare the ampli- 
tude at the rest wavelength of the emission line (designated 
L) to an equivalent width analog: the integral across the 
line's width minus the continuum ratio, divided by L (des- 
ignated W). This is illustrated in figure [TOl This parame- 
terization is useful because it is sensitive to the central am- 
plitude, the integrated strength, and the symmetry of the 
feature, which together are effective in identifying shapes 
similar to the observed microlensing signature. 

Figures [TT] and [T2] show L versus W for a selection of 
the parameter space explored in the simulations. The error 
bars show the location of the data in this parameterization 
and the 2a uncertainty. 

Using this parameterization, we can calculate the p- 
value for the observed microlensing signature with respect 
to each model. The probability, P, of each model reproduc- 
ing the observed microlensing signature within 2a of L and 
W is the ratio of the number of simulated points within 
these error bars to the total number of simulations. The p- 
value is then the integrated probability distribution below 
P, which is equal to the number of simulation points with a 
surrounding point density within a 2cr box equal to or lower 
than the point density of the observed signature, divided by 
the total number of simulation points. Given these p-values. 
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Figure 8. Sample continuum-subtracted line ratios from simulations of outflow model (solid lines), with parameters ro = 1 x 10^^ cm, 
orientation 70° , jet/ disk opening angle 0° /0° . Typical microlensing signatures are asymetric, and the observed hump (dashed lines) is 
rarely reproduced. 
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wavelength (A) 

Figure 9. Sample continuum-subtracted line ratios from simulations of orbital model (solid lines), with parameters ro = 3 x 10^ M©, 
orientation 70°, jet/disk opening angle 0°/0°. The observed microlensing signature (dashed lines) is similar to the simulations in many 
cases. 
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1880 1900 1920 1940 

Rest Wavelength (Angstroms) 

Figure 10. Parameterization of tlie stiape of tlie observed CIII] microlensing signature: we define W , an equivalent width analog for the 
ratio of continuum-subtracted line ratio for images B/A: the integrated ratio minus the continuum ratio {green area) divided by L, the 
central- wavelength amplitude of the ratio {red arrow). 



we can reject models and parameter combinations with a 
quantifyable significance level. 

Figure [T51 shows the p- values for the orbital model as a 
function of orientation, for a range of disk and jet opening 
angles and black hole masses/launching radii. Fig. 1 131 (right) 
shows the average probabilities over all orientations as a 
function of black hole mass/launching radius. 

Figure [14] shows the orientation-averaged p- values for 
both the orbital and the wind models, for all disk and jet 
opening angles and a selection of black hole masses and wind 
launching radii. 

Assuming the orbital model, the p-value for the ob- 
served microlensing signature is inversely proportional to 
black hole mass. We can place a firm upper limit on black 
hole mass, assuming this orbital model. Models with Mbh ^ 
l.SxlO^M© are rejected with >99% confidence. Larger black 
hole masses are ruled out for all orientations and disk/jet 
opening angles. 

The most probable disk/jet configuration for all orien- 
tations and black hole masses is one with a large opening 
angle for both disk and jet of 40° . Interestingly, cases where 
only one of the disk or the jet have a large opening angle 
are ruled out with >99% probability. 

For the pure outflow model, very few parameter sets are 
capable of reproducing the observed microlensing signature. 
All are ruled out with high (> 99%) confidence with the ex- 
ception of models with very small starting radii for the wind 
(10^^ cm; probably unphysical as it is small er than the lower 
estima tes of the accretion disk size, e.g. lEigenbrod et al] 
l2008bl ) and a tightly constrained orientation range (~ 70°), 
which are excluded at only the 90% level. This tells us that 
it is unlikely that a pure accelerating outfiow produced the 
observed microlensing signature. 

Changing the mean microlensing mass marginally af- 
fects the p-values, but the general results are the same. If 
mean masses are smaller, differential microlensing becomes 
less likely but increases in strength, and vice versa for larger 
mean masses. It becomes no easier to produce the observed 
microlensing signature with a pure outfiow model. 



7 DISCUSSION 

Reverberation mapping studies have now provided evidence 
for inflowing, outflowing, Keplerian and virialised motions in 



quasar BELRs (jPennev et aZ.ll2009l : [Sentz et al\ ^200^. Us- 
ing a wide range of observational evidence. iGask cU (200i) 
has argued in favour of infall velocities which are smaller 
than the Keplerian velocities. Flo w-dominated BELR mod - 
els, f or example t hose outlined bv lBlandford fc Reed (|l99ll ) 
and lElvisI (|2000l ) are especially attractive because it is 
known that high velocity nuclear flows are very com- 
mon, if not ubiquitous, in quasars. Indeed, the fraction of 
quas ars in SPSS showing some CIV broad absorption is 
25% (iTrump alf20 06). T he line - driven photoionised wind 
modelled by [Murray et al.\ (| 19951 ') produced line profiles in 
good agreement with observation. In addition, there are now 
strong physical arguments for circular motions (gravity), 
and also the first observational evidence from polarisation of 
the H/3 line (Young et al. 2007). Thus there is still no con- 
sensus on even the simple question of the direction of the 
gas flow in the broad line region. 

Our analysis shows that a pure Keplerian model for 
the BELR is able to reproduce the observed microlensing 
signature. In contrast, our pure accelerating outflow model 
cannot reproduce this signature for any physical launching 
radius. This is expected considering the qualitative nature 
of this signature. The low velocity gas exhibits the strongest 
differential microlensing relative to the continuum (which 
is expected to arise from the smallest size scale), while the 
high velocity gas exhibits flux ratios closer to that of the 
continuum. This distribution of velocities is consistent with 
a Keplerian held in which the highest velocities are expected 
closest to the central gravitating body, but inconsistent with 
a purely accelerating outflow in which velocity increases with 
distance. 

The observed differential microlensing signature for the 
BELR is also consistent with that observed blueness of the 
continuum for image A: if smaller size scales are preferential- 
Uy magnifled, then we expect the continuum to appear more 
blue, assuming that this magnification gradient extends to 
scales smaller the outer regions of the accretion disk. 

7.1 Implications for Outflow Models 

The observed differential microlensing signature 
rules out BELR models in which the gas dynam- 
ics at all radii are dominated by an accelerating 
outfiow. This includes winds that are exclusively 
driven by line pressure. Other drivers, such as ther- 
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Figure 11. "Equivalent widtii" of tiie continuum-subtracted line ratio versus the central flux ratio for simulations {black points) and 
data (red error bars). A selection of the explored parameter space is shown: Mbh of 1 x 10^ Mq to 3 x 10^ Mq, orientations of 30° to 
80°, and jet/disk opening angle of 20°/20°. A black hole mass of 1 X 10^ Mq is favoured, while a black hole mass of 3 X 10^ appears 
unlikely. 
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Figure 12. "Equivalent width" of the continuum-subtracted line ratio versus the central flux ratio for simulations {black points) and 
data {red error bars). A selection of the explored parameter space is shown: ro of 1 X 10^^ cm, 1 X 10^'' cm, and 3 X 10^^ cm, orientations of 
70°, 80°, and 90°, and jet/disk opening angle of 20°/20°. For smaller orientations the data point falls much further from the simulations. 



mal expansion of the confi ning gas (Wevmann et al] 



1982h. centrifugal force jBlandford fc Pavii^ 1 19821 : 
Emmering. Blandford fc Shlosman Il992h or radiation 
pressure via dust (|Voit et al\ 199^ ) may produce outflows 
in which high velocity gas dominates at smaller radii. 

Even if the driving force behind the outflow only allows 
for an accelerating wind, such models are still viable so long 
as the gas remains within the gravitationally dominated re- 
gion for a period comparable to or greater than its orbital 
period. 



7.2 A Pure Keplerian Field? 

Until recently, reverberation mapping studies have found 
that BELR vel ocity distributions are con sistent with Ke- 
plerian motion (jPeterson fc Wandelll200ol ). Our results are 
consistent with this, and extend the finding of reverbera- 
tion mapping — itself applicable to low-luminosity AGN — 



into the high-luminosity regime. Gravitationally-dominated 
kinematics are also supported by the fact that the AGN 
black hole masses calculated based on the assumption of a 
virial BELR follow the Mbh-ct* relation of normal elliptical 
galaxies (Ferrarese et al. 2001). 

Recent reverberation mapping results with greatly 
improved velocity resolution have now found evidence 
for both infall- and outflow-dominat ed BELR kinematics 
(|Dennev et aZ.ll2009l : iBentz et fflZ]|2009l ') , indicating that a va- 
riety of dynamical processes are important. Simple randomly 
oriented orbits are unlikely. It is known that the inner region 
of qu asar BELRs can be as small as 10^® cm jKaspi et al\ 
l200(j ). which overlaps the outer regions of accretion disk. 
This greatly restricts the possible orbits in the innermost 
regions of the BELR. 

Given the likely presence of multiple kinematic regimes 
within the BELR, we interpret our results cautiously. Rather 
than implying a spherical Keplerian field, we instead take 
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these results to indicate that the inner regions of the BELR 
are strongly gravitationally dominated, which may include 
orbital motion and/or infall. 



T.3 A Constrained Vortex? 

In our simple models, the configuration that most easily re- 
produces the observed microlensing signature is one in which 
both accretion disk and jet have large (40°) opening angles; 
configurations in which only the disk or jet have a large 
opening angle are ruled out with high confidence. Interest- 
ingly, this configuration of large opening angles is similar to 
the hoUowed-cone geometry of many wind models. Thus, our 
modeling suggests that such a geometry is viable, but only 
if the gas motion within this outflow cone is gravitationally 
dominat ed at small r adii. This is consistent with the vortex 
model o f lElvid l|2000f ). in which line-driven gas outflowing in 
a hollowed cone initially orbits at accretion disk velocities, 
or with the magnetically-constrained outflow similar to that 
proposed by Emmcring, Blandford & Shlosman (1992). 

This interpretatio n is also consi s tent w ith the biconi- 
cal model proposed by lAbaias et al. 1 l|2007i ) to explain the 
recurrent microlensing enhancement of the broad line blue 
wings in SDSS J1004-f4112. 



7.4 Black Hole Mass Limit 

Our black hole mass limit of Mbh < 1.5 x 10^ M© is of 
similar order to the viri al estimate based on the width of the 
CIV line (0.9 x lO'^MfT,: [Morgan et alhoid ). If we calculate 
a virial black hole mass from the upper limit on the CIII] 
BELR from Section E3] and the line's velocity FWHM, we 
instead obtain Mbh < 2 x 10® M©. 

Modeling the differential microlensing signature with a 
BELR comprised of a simple Keplerian field seems to be 
as useful for obtaining black hole mass limits based on the 
virial assumption as the standard method. It has the added 
advantage of including both BELR velocity and size in the 
model, instead of using quasar luminosity as a proxy for 
BELR size (calibrated using low- luminosity quasars). How- 
ever stronger limits are obtainable by placing a more direct 
limit on BELR size from microlensing of this region. 



7.5 Comparison to Other Works 

lEigenbrod et al] (|2008l ) also detect differential microlensing 
within the velocity structure of the BELR in their long- 
running VLT monitoring campaign. The nature of the signa- 
tures are broadly consistent with those we observe; in high- 
magnification phases, the wings of the CIII] line were more 
strongly amplified than the core. 

Pointdcxtcr fc Kochanekl (|2010l ) are able to constrain 
the orientation of the accretion disk of Q2237-I-0305 by com- 
paring the expected light curve produced by a model accre- 
tion disk with the continuum light curve observed over 11 
years of OGLE monitoring. Their finding of cos i > 0.66, 
which is to say a face-on accretion disk within ~ 50deg of 
the line of sight, is consistent with our results, which don't 
strongly constrain orientation except for a few specific pa- 
rameter combinations; in particular, Keplerian models with 
the largest black hole masses (see Fig. I13p . 
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Figure 13. Left: p-valuc as a function of orientation for a range 
of black hole masses {blue: 10® Mq, green: 3x10® Mq, yellow: 
10^ Mq, red: 1.5x10^ Mq), and tlie average p- value over black 
hole masses {black lines). Right: average p-value over all orienta- 
tions as a function of black liole mass. Vertical panels show re- 
sults for different jet/disk opening angles {upper: 0°/0°, middle: 
20°/20°, lower: 40°/40°). 



8 CONCLUSION 

We have presented observations of the gravitationally lensed 
quasar Q2237-I-0305 taken with the Gemini North CMOS 
IFU. Ratios between continuum subtracted image B and 
image A spectra reveal differential microlensing across the 
velocity structure of the CIII] broad emission line. The high 
velocity wings of this line tend towards the fiux ratio of 
the continuum, and the lower velocity core, while still mi- 
crolensed, is closer to the expected flux ratio in the absence 
of microlensing. This implies that the high velocity compo- 
nent is emitted from a region with a size comparable to that 
of the continuum emission region, whereas the low velocity 
component is emitted from a larger region. 

We conducted microlensing simulations using two sim- 
ple models of the broad emission line region: an outflow 
model and an orbital model. The outflow model assumes 
a clumpy wind accellerated by radiation pressure. The or- 
bital model assumes circular Keplerian orbits with random 
orientations. For both models we tested a wide range of pa- 
rameters. These models were used to construct an ensemble 
of simulated B/A fiux ratios as a function of velocity, for 
comparison with the observed fiux ratio spectrum. 

A purely radial outflow was unable to reproduce the ob- 
served differential microlensing signature for any plausible 
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Figure 14. Probability of the observed microlensing signature 
as a function of model parameters. These values are derived by 
integrating the density functions of the simulated central ratio vs. 
equivalent width distributions (Figs. 1121 and lllll up to the den- 
sity found at the location of the observed signature. Left: orbital 
model, with probabilities averaged over all orientations. Right: 
outflow model, with probabilities averaged over 70° and 80° ori- 
entations only. For all other orientations the outflow model pro- 
duced no results close to the observed signature for any launching 
radius. 



launching radius of the wind. Conversely, the orbital model 
was able to reproduce the observed signature for all sim- 
ulated black hole masses Mi,h < 2 X 10^ Mq. Though our 
orbital model is simplistic, we interpret this result as fur- 
ther evidence that the inner regions of the BELR are grav- 
itationally dominated. This is consistent with an outflow 
model where the BELR gas is lifted off the quasar accretion 
disk, and thus retains a high Keplerian velocity. 

The BELR models presented here are not intended to 
accurately describe the physical situation in the quasar. 
They describe only the generic behaviour of a radially out- 
flowing wind, or a collection of orbiting clouds. More so- 
phisticated models, perhaps making use of radiative transfer 
codes such as CLOUDY, may provide additional constraints 
on the quasar central engine. 

We have obtained Gemini IFU observations of nine 
other gravitationally lensed quasars, both double and 
quadruply imaged. These data will allow us to probe a range 
of quasar orientations, emission lines, black hole masses, and 
therefore BELR emission region scales. The analyses of these 
data are forthcoming. 
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